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Abstract 

Nonrenormalizable quantum field theories require counterterms; 
and based on the hard-core interpretation of such interactions, it is 
initially argued, contrary to the standard view, that counterterms 
suggested by renormalized perturbation theory are in fact inappro- 
priate for this purpose. Guided by the potential underlying causes of 
triviality of such models, as obtained by alternative analyses, we fo- 
cus attention on the ground-state distribution function, and suggest a 
formulation of such distributions that exhibits nontriviality from the 
start. Primary discussion is focused on self-interacting scalar fields. 
Conditions for bounds on general correlation functions are derived, 
and there is some discussion of the issues involved with the contin- 
uum limit. 
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1 Introduction 



Nonrenormalizable quantum field theories have long been the "black sheep of 
the family".^ While supcr-rcnormalizable and strictly renormalizable quan- 
tum field theories (especially asymptotically-free theories) have been well 
served by renormalized perturbation theory, this has not been the case for 
the usual nonrenormalizable family members. For them - so the pertur- 
bation story goes - an infinite set of distinct counter-terms are necessary, 
requiring thereby an infinite number of experiments to establish the needed 
coefiicients of the counterterms, an endeavor which nullifies any predictive 
power the theory may have had. A key concept in the preceding sentence re- 
lates to the perturbative counterterms, and we shall stress the shortcomings 
of this concept below. 

Nonrenormalizable models may also be approached nonperturbatively. 
Here, we have in mind the analysis of v?^ models for spacetime dimension 
n > 5 made by Aizenman [1] and Frohlich [2]. In their approach, they for- 
mulated the models as Euclidean-space lattice theories with arbitrary coeffi- 
cients for the mass term and the coupling constant, both taken as functions 
of the cutoff parameter, namely, the lattice spacing a; of course, the coupling 
constant was required to be nonnegative. They were rigorously able to show, 
independent of the choice of the indicated parameters, that the continuum 
limit (as combined with the infinite volume limit) of the theory led to the 
result of a (possibly generalized) free field. In other words, the manifestly 
non-Gaussian nature of the lattice field distribution has, in the continuum 
limit, become a strictly Gaussian distribution. Since Monte Carlo calcula- 
tions suggest that triviality holds for iff as well, this overall story has given 
rise to the widely held view of the "triviality of the </?^ theories" . 

In the author's judgement, the infinite number of distinct counter- terms 
as well as the triviality result, while correct conclusions given their respective 
underlying assumptions, are both unacceptable as physically correct quan- 
tum statements for the problem at hand. This conclusion is based on the 
following assertions: (1) For a nonrenormalizable theory, results based on 
regularization and the introduction of counterterms given by a perturbation 
analysis al)out the free theory are unacceptable for the simple reason, as 

^Common meaning: "A worthless or disgraced member of the family" . See, for example, 
http://www.phrases.org.uk/meanings/66250.html. 
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we shall argue, that the interacting nonrenormalizable theories are not even 
continuously connected to the free theory as the coupling constant is reduced 
to zero; (2) The classical (pf^, n > 5, models are manifestly nontrivial theo- 
ries exhibiting nonvanishing scattering, etc. If the quantum formulation of 
such models is taken to be the trivial one, then it follows that the classical 
limit of that quantum theory leads to a trivial classical theory which is not 
equivalent to the original classical theory one started with. Such a dilemma 
can only mean that the quantum formulation leading to triviality cannot be 
the physically correct quantum formulation for nonrenormalizable </?^, n > 5, 
models. 

We claim there is an alternative formulation of the quantum theory for 
such models that can explain the unacceptable results of infinitely many 
counterterms as well as triviality, and which furthermore offers alternative 
calculational possibilities that we believe may lead to satisfactory results. 
The purpose of the present paper is to lay out the theoretical arguments 
supporting the point of view under present consideration. We start by reex- 
amining the limiting Gaussian behavior of the ip'^, n > 5, models. 

In probabihty theory, the Central Limit Theorem (CLT) leads to Gaus- 
sian distributions for a very large class of independent, identically distributed 
random variables. By concentrating attention on random variables that are 
independent and identically distributed, it is of course possible to demon- 
strate the CLT quite generally. However, it stands to reason that there exist 
distributions of random variables that are neither identical nor independent, 
but which are close to such behavior in some unspecified manner, yet the 
combination nevertheless tends to a Gaussian distribution as the number of 
variables increases without limit. General theorems for such cases arc natu- 
rally most unlikely since a Gaussian or non-Gaussian limit depends explicitly 
on the particular details at hand. Nevertheless, based on the evidence, it is 
plausible to regard </?^, n > 5, models as satisfying the criteria to lie within 
the basin of attraction of the CLT. Our further development will be partially 
infiuenccd by this viewpoint. 

In Sec. 2 we outline the argument why nonrenormalizable interactions 
act as discontinuous perturbations, and as such, counterterms suggested by 
regularized and renormalized perturbation theory are unsuitable. Attention 
shifts from the lattice action and the field distribution it engenders in Sec. 3, 
and instead focusses on the sharp-time, ground-state distribution. Arguing 
from the CLT, we propose that this distribution should be a generalized 
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Poisson distribution, a subject discussed at length in Sec. 4. A generalization 
to a linear superposition of Poisson distributions is the subject of Sec. 5. It is 
shown in Sec. 6 that general spacetime correlation functions can be expressed 
as, or are at least bounded by, properties of the Poisson distribution, while 
Sec. 7 offers a useful reformulation of some of the expressions developed in 
Sec. 6 as well as a brief discussion of the continuum limit. Section 8 is 
devoted to a conclusion, and some details about specific properties of a given 
generalized Poisson distribution are discussed in the Appendix. 

It is not without interest that in an archived (but unpublished) paper [3], 
the author already made the proposal that the ground-state distribution for 
nonrenormalizable scalar models should be a generalized Poisson distribution. 
Unfortunately, how that proposal was used in that earlier paper to suggest 
an auxiliary potential and thus a lattice action was incorrect. Hopefully that 
situation has been advanced in the present paper. 

2 Inadequacy of Regularized, Renormalized 
Perturbation Theory 

2.1 The hard-core picture 

Chapter 8 in [4] describes in some detail the hard-core picture of nonrenor- 
malizable interactions, and for a detailed account of that approach we can do 
no better than to encourage the reader to examine that presentation, which, 
importantly, is for all intents and purposes an independently readable chap- 
ter. However, for the sake of the reader who does not need a detailed account 
of the hard-core picture of nonrenormalizability, we offer the following brief 
synopsis. 

In general terms, let the theory in question be described by a formal 
(Euclidean-space) functional integral commonly called a generating func- 
tional, which for the "free" theory is given by 



where /i is a smooth "source" function, and Wo(0) > denotes the free 
field (Euclidean) action, an expression which is often quadratic in the fields 
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but need not be so. Likewise we have the generating functional for the 
"interacting" theory 



Sx{h} = AfJ e/'^'^'^"--^°(^)-^^('^) V<f> , 



where the couphng constant A > 0, and V{(j)) > denotes the (Euchdean) 
interaction action. In each case, the formal normalization is chosen so that 



for all h (say C^) ; in addition, we say that y is a discontinuous perturbation 
of Wo if, instead. 



for all h, where 5'o{/i} is referred to as the generating functional for the 
"pseudofree" theory. 

At first glance, it seems self evident from the given expressions for the 
generating functionals that the limit of Sx{h} as A — > must be So{h}. 
Moreover, the unquestioned belief that this must be so would appear to be 
a tacit assumption behind any perturbative analysis of Sx{h} in a series 
expansion in A about A = 0. However formal (i.e., nonconvergcnt) that series 
may be, the introduction and use of such a series expansion presumes that 
y is a continuous perturbation of Wq. If V is not a continuous perturbation, 
i.e., if V is in fact a discontinuous perturbation of Wq, then there is no 
way in which one could legitimately claim that ^^Sx{h} has a perturbation 
series expansion about So{h}" . That statement is as likely to be true as the 
statement that properties of the (discontinuous!) function 



can be determined by a power series expansion about A = 0. 

Qualitatively speaking, a continuous perturbation is one for which 




lim Sx{h} = So{h} , 

A — >■© 



hm Sx{h} = S',{h} ^ So{h} 

A — *U 



F(A) = e^, A > , 

F(0) = TT 
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Similarly, a discontinuous perturbation is one for which 

Wo((l)) < oo ^ < oo , 

or in other words, there are fields (comprising a set of nonzero measure) for 
which 

Wo{(l)) < oo , = oo . 

For a discontinuous perturbation, the contributions such fields would have 
made if the interaction term V had never been present are now missing 
altogether because those fields are effectively projected out by the factor 
e~^^, for any A > 0, however small. If we introduce 

= 1 , if W((t)) < oo , < oo , 

= , if W((l)) < oo , = oo , 

then we see that we can freely redefine Sx{h} as 

Now, as A —> 0, we can readily see that we obtain 

S'^{h] = M'o j e/^*'^"^-^°W X(0) V4) . 

In brief, the interaction term has acted partially as a hard core projecting 
out - thanks to X{(1)) - certain contributions that would have been included 
for the free generating functional! 

To show the relevance of this discussion to 0^ models, we need only recall 
for n < 4 (the renormalizable cases) that 

{/0^(x)d"x}V2 < (4/3)/{[V0(x)]2 + m2(/)2(x)}d"x 

for any w? > 0, while for n > 5 (the nonrenormalizable cases) there are 
singular fields, such as 

e n n 
^ x\P '4^2 
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for which 



/{[V0,(x)]2 + m'0^(x)}d"x < oo , 

while 

(tx ^ oo , 

fulfiUing the requirement for the interaction to act as a hard core. 

A derivation of the (4/3) bound given above when n < 4 is part of the 
discussion in Chapter 8 of [4]. The factor (4/3) also provides a bound for 
the more general models when p < 2n/{n — 2) (the renormalizablc cases), 
while there is no finite bound when p > 2n/(n — 2) (the nonrenormalizable 
cases). When the gravitational action is divided into quadratic and non- 
quadratic terms, the latter term has the properties to act as a discontinuous 
potential consistent with the well-known nonrenormalizability of gravity. In 
fact, as also detailed in Chapter 8 of [4], discontinuous perturbations of the 
harmonic oscillator arise even in conventional quantum mechanics. If hard- 
core potentials exist in quantum mechanics, it should be no surprise that 
they arise in quantum field theory! 

We are about to leave the general description of hard-core interactions 
hoping that the reader has accepted that hard-core interactions may indeed 
exist and that nonrenormalizable quantum field models are likely candidates. 
Harder to accept, but nevertheless a direct consequence of the existence of 
hard-core interactions, is the concept that in such cases counterterms sug- 
gested by a (regularized) perturbativc power series expansion about the orig- 
inal free theory are not relevant - indeed, harsh as it may seem, we can even 
say they are completely irrelevant! In fact, such terms are entirely mislead- 
ing because regularized, perturbative counterterms are designed to preserve a 
continuity of the interacting theory with the original free theory, a continuity 
which in fact never existed! 

This statement about perturbation theory is no doubt difficult to accept 
for some readers because it certainly goes against - even contradicts - the 
"general wisdom" of how nonrenormalizable models are normally presented. 
If the reader's prejudice in favor of perturbation theory and its unquestioned, 
universal applicability prevent you from accepting the statement that in gen- 
uine hard-core cases perturbation theory about the original free theory is ir- 
relevant, then this article is not for you. By all means, stop reading here and 
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now! However, if you can accept the fact that hard-core interactions have 
the power to change almost everything regarding perturbation theory about 
the free theory, then we invite you to read on. Indeed, we can offer you the 
comforting possibihty that the interacting theory may indeed have a mean- 
ingful perturbation expansion - not about the free theory - but about the 
pseudofree theory, which after all is the theory to which it is continuously 
connected as A — > 0! As an elementary illustration of that very thought, 
let us revisit the discontinuous function F{X) introduced earlier. This time, 
however, we introduce an analogue of the pseudofree theory by letting 

F'(A) = F(A) = e\ A > , 
F'(0) = limF(A) = 1 . 

Clearly - and unlike the original function F(A) - the new function F'(A) 
does have a pcrturbative power series representation about the "pseudofree" 
value, -F'(O) = 1, namely, 

n=0 

which is valid in the present case for all A > 0. More importantly, highly 
specialized, soluble nonrenormalizable models (see Chapters 9 & 10 in [4]) 
possess solutions that do admit entirely reasonable perturbative treatments of 
their interactions about the appropriate pseudofree solution. We believe valid 
perturbative analyses about suitable pseudofree solutions hold widely, and we 
will make the assumption in this paper that such perturbation formulations 
exist. 

Accepting the concept that the nonlinear interaction for <^^, n > 5, mod- 
els acts partially like a hard core leads immediately to the conclusion that 
some form of counterterm is needed, at the very least, to represent the ir- 
removable effects of the hard core. The triviality argument also supports 
the need for counterterms since mass, coupling constant, and field-strength 
renormalization alone are insufficient to escape triviality. Our previous dis- 
cussion has been necessary to open the way to consider counterterms other 
than those suggested by renormalized perturbation theory. Our task is to 
find suitable counterterms that will avoid triviality - a highly nontrivial task 
to be sure! 
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3 Focus on the Sharp-Time, Ground-State 
Distribution 



3.1 Triviality also holds at sharp time 

A Euclidean spacetime distribution that exhibits triviahty is necessarily Gauss- 
ian, and has vanishing truncated correlation functions beyond second order. 
Such correlation functions a fortiori vanish at equal times, a fact that im- 
plies that the associated ground-state distribution for the limiting Hamilto- 
nian has a Gaussian distribution itself. Stated otherwise, the lattice-space 
ground-state distribution - the square of the lattice-space ground-state wave 
function itself - which because of the nonlinear interaction is manifestly non- 
Gaussian, has, in the continuum limit, become Gaussian. Just as the full-time 
distribution has become Gaussian because it lies in the appropriate basin of 
attraction for the CLT, it is equally plausible to conclude that the sharp-time 
field distribution for the nonlinear theory also lies in the associated basin of 
attraction of the CLT. Conversely, if we can ensure that the ground-state 
distribution remains non-Gaussian, then the non-Gaussian character of the 
full-time distribution will be assured. This remark leads us to focus some 
attention on the putative ground state. 

3.2 Lattice space version of the ground-state wave func- 



We let the even, positive function ^{(p) [— ^'(— 0)], denote the ground-state 
wave function on the Euclidean-space lattice, and we adjust the potential 
by a constant so that this ground state has zero energy. Consequently, the 
nonnegative Hamiltonian operator for the lattice is given by 



tion 




where the potential is given by 
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Let us assume that the lattice has lattice spacing a, is hypcrcubic with pe- 
riodic boundary conditions in all coordinate-axis directions, and the number 
of lattice sites k — {ko, ki, k2, ■ ■ ■ ,ks), kj e Z, s — n — 1, in each axial di- 
rection is L. Here the primed sum E'^ signifies a sum over spatial values of 

I.G.J ^2; • • • ; 

some specific, but generally implicit, fixed, tempo- 
ral value of k, i.e., kc,. The argument stands for the collection of lattice 
fields {(pk} at an equal lattice time. Thus if /c = {ko, ki, k2, ■ ■ ■ , kg), we se- 
lect ko as the future time direction (under a possible Wick rotation) , and so 
4> = {(l>k '■ all k with ko fixed}. 

The probability distribution constructed as the square of the ground- 
state wave function may alternatively be characterized by its characteristic 
function 

C{h) = j e'^'"^"^^"' ■^{(t>fXl'j^d(t)k] 

the symbol IT'^ denotes a product over spatial values of k holding the temporal 
value fixed. We define the continuum limit of such expressions to mean: (i) 
the lattice spacing a — > as well as (ii) the spatial volume {Lay — > oo 
along with the spacetime volume {LaY ^ oo in an appropriate manner and 
combination. As noted above, if we choose the naive lattice form (including 
arbitrary cutoff dependent parameters) for the classical action of a > 5, 

model, the ground-state distribution becomes Gaussian, or equivalently that 

C{h) Eoaussianih) = ex.p[- h(x) U (x - y)h(y) d^xd^] , 

for some covariance function U{x — y). Our primary goal is to avoid a Gaus- 
sian limiting behavior in the continuum limit for the ground-state distribu- 
tion. 



3.3 Some partial matrix elements 

The lattice-space Hamiltonian discussed above admits an alternative and 
useful form. Let us introduce the differential operator 

Afe = - a'^ — + a-' — - — ^ , 
d(pk W(<p) dcpk 
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and its adjoint operator 

for each A; in a given spatial lattice slice. It is straightforward to verify that 
the lattice Hamiltonian is also given by 

Next consider states of the form 

and note that 

Consequently, it follows that 

J ^h'{<PTn^h{<P) Kd<Pk = h'.hk C{h - h') . 
Additionally, we observe that 

4(4-^,) = — 
the canonical momentum conjugate to the field (pi. It follows that 

*^,(0)*7ri*ft(0) n',#fc = m + hi) C{h - h') . 



Equivalently, if we use the connection 

d 

m = -ia = i[H, (pi] 
odi 

relating the canonical momentum with the field and Hamiltonian, then it 
follows that 



/ 



\{h'i + hi)C{h-h') , 
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as before. (Remark: There are now THREE different uses of the prime: one 
refers to the pseudofree theory; another refers to properties on a fixed time 
shce of the lattice; the last usage refers to labels of Hilbert space functions. 
The context suffices to tell which usage is meant.) 

Observe that this calculation has yielded matrix elements of both the lat- 
tice Hamiltonian Ti and the lattice momentum tt; at site /. It is highly proba- 
ble that these matrix elements actually determine the operators in question, 
albeit in an indirect manner. 

3.4 Features of the continuum theory 

The analysis given above for the lattice has a natural analog in the continuum 
theory [5]. The expressions given below should indeed follow directly as 

the continuum limit of the lattice expressions, despite the fact that wc have 
changed notation to emphasize the fact that we are dealing with the putative 
continuum theory. 

The continuum limit of the ground state may be formally identified with 
the abstract unit vector |0), a member of the abstract Hilbert space ^. 
Additionally, we can introduce the unit vectors 

\h) = Qi'J'^(^)4>{x)d''x |Q^^ ^ 

where s = n — 1, 4>{x) denotes the sharp-time (e.g., at t = 0), formally 
self-adjoint field operator, and h{x) denotes a smooth, real test function. 
Additionally, the vector \h) is strongly continuous in a suitable topology [6]. 
The overlap of two such vectors, 

E{h-h') = {h'\h) , 

defines the important expectation functional E{h). As a continuous function 
of positive type, E(h) serves as a reproducing kernel for a reproducing kernel 
Hilbert space that provides a representation of by continuous functionals. 
Furthermore, if we define the local operator 7r(x) as the canonical conjugate 
of the field operator 0(x), then, under mild conditions, it follows [5] that 

{h'\Ti{x)\h) = l[h'{x) + h{x)] E{h - h') . 

Since we may assume that the vectors \h) span the Hilbert space, it is plau- 
sible that the given matrix elements of ti{x) determine the local self-adjoint 
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operator uniquely. Moreover, we can also assert [5] that 

{h'\n\h) = \{h',h) E{h-h') , 

where Ti. denotes the nonnegative, self-adjoint Hamiltonian operator for which 
7i|0) =0, provided we assume that 

7t{x) — i[H, (f){x)] . 

Note well, that these relations which hold in the continuum are simply 
natural continuum analogs of the corresponding lattice-space expressions de- 
rived above, and confirm that a smooth continuum limit of the function C(h) 
to E{h) is a primary ingredient in developing a satisfactory continuum theory. 

4 The Role of Generalized Poisson 
Distributions 

In probability theory, alternative distributions of the sum of an ever- increasing 
number of independent, identically distributed random variables that com- 
pete successfully in the struggle for convergence, such as in a continuum 
limit, generally lead to Poisson distributions, and in doing so they give rise 
to non-Gaussian final distributions. Significantly, a sequence of lattice-space 
ground state distributions with a limiting behavior that leads to a Poisson 
distribution necessarily has a qualitatively different dependence on the indi- 
vidual variables and various parameters than does a sequence of lattice-space 
ground state distributions destined to end up as a Gaussian distribution. 
Just as for the Gaussian limit, it is plausible that there are distributions for 
alternative sets of random variables, which are neither independent nor iden- 
tically distributed, yet lie in suitable basins of attraction so that they have 
limiting Poisson distributions. If we ensure that the terms in the lattice ac- 
tion along with appropriate counterterms conform with the needed different 
dependence on the individual variables and parameters, then we can ensure 
that the continuum limit will avoid being Gaussian. 

As stressed previously, the traditional free theory is Gaussian and there- 
fore corresponds to an infinitely divisible distribution (defined below) . There 
is some logic in the fact that a free theory is infinitely divisible in the sense 
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that the abihty that each field variable may be expressed as a sum of ar- 
bitrarily many independent, identically distributed variables is, let us say, 
one hallmark of being "free" . Likewise, when discussing a pseudofree theory 
in which all interactions have been reduced to zero, it is not unrealistic to 
believe that the hallmark of infinite divisibility also applies to the resultant 
pseudofree distribution. While arguments of any weight are hard to advance, 
we shall nevertheless make the assumption that it is possible that this is 
the case and, as a consequence, we shall seek pseudofree solutions among 
the family of infinitely divisible distributions. Even so, we shall also hedge 
our bets by suggesting that the class of ground-state distributions which are 
given by rather general linear combinations of infinitely divisible distribu- 
tions (and are therefore, in general, not infinitely divisible themselves) may 
be the proper place to find pseudofree theories. 

Before studying particular examples of infinitely divisible distributions, 
however, it is appropriate to give a general outline of how such distributions 
may be characterized. 

Successful limiting behavior of suitable sets of countable numbers of in- 
dependent, identically distributed random variables may be formalized in a 
natural way in terms of their characteristic functions. In particular, the ap- 
propriate Gaussian and Poisson distributions share the property of infinite 
divisibility, which may be described as follows: If C{h) denotes an infinitely 
divisible characteristic function for an appropriate distribution, then it fol- 
lows that the J^^ root of that function, i.e., Cihy/-^, J e {2,3,4,...}, is 
also a characteristic function. This property follows because we can always 
linearly decompose the set of random variables {0^} into J components com- 
posed of independent, identically distributed random variables for any pos- 
itive integer J. Furthermore, it is always the case that any integer power 
X of a characteristic function is again a characteristic function, i.e., if C{h) 
is a characteristic function, then C{h)^, K e {2,3,4,...} is also a char- 
acteristic function. Therefore, an infinitely divisible characteristic function 
C{h) retains the property of being a characteristic function if we raise it to 
an arbitrary rational power such as C{h)^^'^. Finally, since characteristic 
functions are always continuous functions of their arguments, we can take a 
limit as the rational ratios converge to an arbitrary, positive real, K/ J ^ r. 
This limit converges, and the limit C{hY is again a characteristic function 
for an arbitrary real number r > 0. This feature is a defining as well as a 
useful property of infinitely divisible characteristic functions and thereby of 
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infinitely divisible distributions. 

The only categories of infinitely divisible distributions are either Gaussian 
or Poisson, or a combination of the two. And since a Gaussian distribution 
can be obtained as the limit of a sequence of Poisson distributions (but not 
vice versa!) it is appropriate to say that all infinitely divisible distributions 
are Poisson distributions or limits thereof. 

In what follows we shall assume that the ground-state distribution func- 
tion has the form of a fairly straightforward, basic Poisson distribution. We 
conjecture that some of the basic Poisson distributions on which we focus 
could possibly serve as appropriate pseudofree scalar field models. 

In Sec. 5, ground-state distributions that arc given by linear superposi- 
tions of basic Poisson distributions are treated as possibly interesting gener- 
alizations of the case of simple Poisson distributions. 

In seeking models within the class of strictly infinitely divisible distribu- 
tions, or within the class of linear superpositions of such distributions, it is 
important to keep in mind that some specific model, such as n > 5, may 
not be among those contained in the class under consideration. Indeed, one 
should not look for any specific model, but, at the present stage, accept any 
reasonable model that might be present. This view is necessitated because 
we are faced with a highly nontrivial inverse problem; namely, we can access 
a limited amount of information about the model under consideration, i.e., 
some features of the ground-state distribution, and from that we would like 
to determine the lattice action that gave rise to those properties. In sum- 
mary, we ask the reader to keep an open mind as to just what class of models 
of interest - if any - can be described by the procedures we have in mind. 

4.1 Basic Poisson distributions 

To describe a Poisson distribution for scalar field models in terms of the 
ground-state distribution ^(0)^ is difficult because it generally does not in- 
volve known functions. On the other hand, the characteristic function for a 
general Poisson distribution can be described in a relatively simple form. In 
particular, the characteristic function of an even Poisson distribution has the 
form given by 

J cos{E'khk(l>ka') n'f^dcPk = exp{-/[l - cos{E'khk(l>ka')]p{(f), a) U'^dcPk} , 
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where the weight function p(0, a) > 0. 

To ensure absolute continuity of the ground-state distribution, it is nec- 
essary and sufficient that 



Jp{(j),a) Il'f^d(j)k = oo , 



a criterion that estabhshes our choice of Poisson distributions as so-called 
generalized Poisson distributions [7] . And to ensure proper meaning of the 
right-hand side of the equation for the generalized Poisson characteristic 
function above, it is necessary that the exponent - frequently referred to as 
the "second characteristic" [8] - satisfies 



for all suitable {hk}- By symmetry, the weight function p{—(j),a) = a), 
and we assume that all even "moments" exist in the sense that 



for all q > 1, where /c™ = (/cq, k"\ klj^, . . . /c™), and the time ko is the same for 
all A;™, 1 < m < 2q. By assumption all odd "moments" vanish. Note that 
such expressions arc not moments since the weight function p{(j), a) is not - 
and cannot be - normalized. To emphasize that distinction we shall generally 
refer to such expressions as noments (rhymes with moments, and is derived 
from not moments, a phrase that determines the meaning of noments). 

In order that p(0, a) have the indicated properties, it must be singular 
near = in a suitable manner. For our purposes we shall focus on a) 
of the form 



Let us explain the separate terms that enter this expression. 

The numerator contains a factor R{a) to be fixed later and an unspecified 
function f/(0, a), which is limited at present, let us say, by the requirements 
that U (0, a) = and U (0, a) > cT,'j.(j)la'^ , for some c > 0, and which is defined 
for field variables at a fixed time ko; at the present moment, the only purpose 
of U {(f), a) is to ensure that all noments of p{(f), a) are finite. Below, we shall 
further specialize U {(f), a) to the quadratic form 



< /[I - cos(S'^/ifc0fca*)]p(0, a) Tl'^d(f)k < oo 





U{(f),a) = i:'^i(f)kAk,i(f)ia 
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for suitable (but unspecified) matrices Ak^i for closer examination as candi- 
dates for pseudofree models; we shall even take the liberty of calling models 
with this form "pseudofree models" . As will become clear, one argument for 
this designation is the fact that in its final form, this particular version of 
the weight function p(0, a) has the same number of free, dimensional param- 
eters as a traditional free theory. Additional arguments in favor of such a 
designation will arise later. 

As shown in the Appendix, the exponent 7 in the denominator of the 
expression for a) satisfies the bounds 



1 11 

- < 7 < - H 

2 - ' 2 A^' 



where N' = = L""^ denotes the number of lattice sites in a spatial 
slice. Since we eventually need to take A"' — > 00, we accommodate that limit 
already by choosing 

1 

the only exception being in the Appendix where we establishing the bounds 
on 7 noted above. Since the parameters {Pk} will turn out to be dimension- 
less, then, when 7 = 1/2, it follows that R{a) is dimensionless. Consequently, 
all the dimensional parameters in p{(f),a) reside in U{(f),a). 

The factors (5k [= P-k] are all nonnegative, i.e., ^ 0, and satisfy the 
condition that 

KPk = l. 

Various choices of the set {Pk} are possible, but we shall concentrate on 
choosing a fixed, A^'-independent number of (3k terms as nonzero. In partic- 
ular, we specialize to the specific choice of 2s -|- 1 nonvanishing terms given 

by 

Pk ^ Jk= ■^^^-^^k,{kUk„„} ■ 

This notation means that an equal weight of l/(2s -|- 1) is given to the 2s -|- 1 
points in the set composed of k and its 2s nearest neighbors in the spa- 
tial sense only. Specifically, we define = I /{2s + 1) for the points k = 
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{ko, ki, k2,...,ks), k= {ko, /ci ± 1, /c2, . . . , kg), k = (fco, ki,k2±l, . . . , ks),. . . , 
k = {ko, ki, k2, ■ ■ ■ , kg ± I), and = for all other points in that spatial 
slice. ^ The suitability of this choice is addressed in the Appendix. For nota- 
tional purposes we shall refer to the set of lattice sites {k U knn} as the set 
The name /c^'^^ is chosen because for the graphically simple case for 
which s = 2, the set of nonvanishing {Jk} points arc arranged in the form 
of a "+" sign; indeed, see Fig. 1 in the Appendix in this regard. (Remark: 
Of course, one could give an unequal, positive weighting to each of the /c^'"^ 
points in contrast to the set of equal weights we have chosen. While this 
choice would lead to a different sequence of results for finite lattice spacings, 
any difference due to unequal weighting would disappear in the continuum 
limit. Such an argument applies to any set {Pk} that contracts to a single 
point in the continuum limit.) 

Clearly, not all ground-state distributions ^(0)^ correspond to infinitely 
divisible distributions. However, for every acceptable choice of U{(f),a), the 
prescription that defines C{h) automatically defines - albeit indirectly - a 
ground-state distribution \1'(0)^ that is infinitely divisible. 

4.2 What is the pseudofree ground-state distribution? 

Given a characteristic function for a generalized Poisson measure, the asso- 
ciated distribution is given be an inverse Fourier transform, namely. 



X exp{ — J[l — cos(E'fc/ifcMfca*)] p{u, a) 11'f.duk} H'kdhk ■ 

In general, this integral cannot be evaluated analytically; however, there 
is a partial result which is illuminating. It follows from the property of 
infinite divisibility that C{hY^^ is again a characteristic function for all 
M G {2, 3,4,.. .}, and we can put that property to use for us now. 
Consider the expression 



^Originally [3], it was believed that it was necessary that all /?fe > was required. We 
thank Erik Deumens for raising the question of whether finitely many positive jSk terms 
might work just as well as requiring that all /3fe > 0. Happily, that has turned out to be 
the case. 
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X 



exp{— (1/M) /[I — cos(E'^/ifctifca*)] p{u^ a) Ul^^duk] H'kdhk , 



which yields the distribution belonging to the characteristic function C{hy^^. 

For very large M it stands to reason that the second characteristic is small, 
which allows us to approximate the foregoing integral up to order 0(1/M^) 
as follows: 



s s N' 



1/M 
X 



V27r 



cos(Efc0jk/ifca*) 

1 - (1/M) /[I - cos{E'f^hkUka')] p{u, a) U'^duk] U'^dhk , 

Now we make use of the specifically assumed form for a) given by 

R{a) exp[ — ?7(0, a)] 



p((/),a) 



and remark, since we are working to order 0{1/M) that we can replace this 
weight function by (say) 

_ R{a) exp[-U{(j),a)] 
^'/^^'^''^^ - n'jSi J.-,0? + F(M)]V2 ' 

where the additional term in the denominator, F{M), satisfies the bounds 

< F(M) < 1 . 

With this replacement in place, we now claim, with the same accuracy as 
before, that 



s X TV' 



2n 



cos(Efc0fe/ijka*) 



X 



1 - (1/M) /[I - cos(E^/ifc'Ujka'')] Pi/m{u, a) Il',^duk U-'kdhk . 



We are permitted to make this small change of introducing F{M) because our 
normal integrand for the second characteristic is "protected" from diverging 
for very small values of by the factor 



1 - cos(E^/ifc0fca'') 
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in the numerator. However, - and unlike a) - the weight function 
Pi/m(0, is integrable, and we now adjust F{M) so that the integral is 
M itself; i.e., we now assume that F{M) is chosen so that 

J Pi/m(0, a) W^d(j)k = M . 

As M — s> oo, it is clear that F{M) — >■ 0; however, any further properties of the 
function F{M) are not of particular interest to us. With the normalization 
of Pi/m(0, fixed, we may recast the expression for ^^^^^^(0) into the form 

*i/M (</>)' = (1^) / cos(E'fe(/)fc/ifca^) 

(1/M) J cos{^'khkUka')pi/M{u,a) n'^diife] H'^d/ife , 



X 



an equation which permits us to identify 

vi> r^^2 ^ J_ Rja) exp[-[/(0,a)] 

correct to order (9(1/M^). 

This latter equation makes clear, for large M, that the ground-state dis- 
tribution for this special case is just the weight function a) with its sin- 
gularity regularized and then rescaled to have integral one. Apart from the 
denominator factor, the ground-state distribution ^'i/m(0)^ is proportional 
to exp[—U{(f), a)], which for the particular case in which U{(f), a) is quadratic 
in the fields is highly suggestive of a free model. Here is further evidence 
supporting our designation of weight functions of the form 

p[(p,a) = — 



as candidates for pseudofree models. (Remark: It is noteworthy that cer- 
tain highly specialized, soluble, nonrenormalizable models (see Chapters 9 
& 10 in [4]) also exhibit ground-state distributions of a qualitatively simi- 
lar form, including certain denominators not unlike our featured expression 
H'jS;jfc_«0^ + F(M)]i/2. In those cases, the analog of our function U (0, a) 
was necessarily quadratic for the associated pseudofree models.) 



20 



While the characteristic function of interest is easily regained as the Mth 
power of the characteristic function C{h) ^^'^^ , it is not nearly so easy to regain 
the ground-state distribution of interest since ^(0)^ is given by the M-fold 
convolution of the distribution '^i/m{4'Y-i followed, ultimately, by the limit 
M ^ oo to eliminate the error made in expanding the exponent. If we let the 
symbol * denote convolution, then the computation just outlined is given by 

= lim *i/m(0)^ * *i/m(0)^ * • • • * *i/m(0)^ , 

M — ^oo 

where the multiple convolution involves M factors. 

After multiple convolutions, what happens to the original form of the 
expression for the ground-state distribution ^i/m(0)^ is hard to predict. If 
U{(f),a) is chosen as a quadratic function, which we have advocated for our 
pscudofrcc models, and if we momentarily assumed that the denominator in 
p{(f), a) were absent, then the original distribution is a pure Gaussian, and 
under repeated convolutions it would still remain Gaussian; this is the usual 
situation for a strictly free model. However, the denominator of a) is ex- 
actly what separates the Gaussian and the generalized Poisson distributions, 
and therefore the denominator factor is of the utmost importance. 



4.3 A glimpse at the potential 

We are unable to offer an analytic form for the ground state ^(0), but we 
are able to offer an analytic form for the ground state ^i/jv^(0) valid for large 
M. Of course, the latter ground state is not the one of ultimate interest, 
but for someone starved for a look at any form of a ground state, it may 
offer a ghmpse into the problem that is unavailable elsewhere. The potential 
associated with the ground state ^i/m(0) is given by the general formula, 
namely 



i/m(0) d<Pl 
which leads to the expression 

1 -sV^' Js-kJt-k(t>i 



Vi/M(0) = -a-J]' 



8 [T^'mJs-mCkl + F{M)\ [KJt-n<l>l + F(M)] 



21 




2 + Vi/m(0) , 



k,t 



where V-\_/m{4>) follows from terms which include either one or two derivatives 
of the unspecified term C/(0, a). If we choose 



it follows that the remainder potential Vi/m{4>) contains a quadratic poten- 
tial term plus a term that arises from one derivative of U and one of the 
denominator. Observe well, that when factors of H are taken into account, 
the quadratic term is 0{h^), the cross term is 0{h}), while the term displayed 
in Vi/m(0) above - coming from two derivatives of the denominator factor - 
is 0{h'^). In the classical limit, therefore, the only term that survives is the 
quadratic term; the other two terms are quantum corrections to the potential 
that arise strictly from the denominator factor present in the weight function 
Pi/m{(P,ci), and they make no contribution in the classical limit. Since the 
denominator carries no dimensional parameters, there is no possible way for 
it to contribute to the classical limit! 

Of course, we must not forget that we are only looking at the potential 
associated with \E'i/m(0) and not that associated with ^(0), the one of real 
interest. There is no simple way to tell how the potential V(0) will appear af- 
ter the multiple convolutions that ^i/m(0)^ must undergo to become ^^(0)^; 
nonetheless, there is no denying the feeling that something like the contri- 
bution displayed in the expression for Vi/m{4>) will survive, if for no other 
reason, than on dimensional grounds alone. However, that is only a guess. 

Understanding repeated convolutions of the distribution ^i/m(0)^ for 
large M is one key to getting a handle on the ground-state distribution 
^(0)^ of ultimate interest. Nevertheless, there may be another way to gain 
some information about ^'(0)^, and that is to proceed numerically. 

4.4 Possible numerical studies 

Whatever expression one chooses for the characteristic function C{h) of the 
ground-state distribution ^'(0)^, it follows that the ground-state distribution 
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itself is given simply by the inverse Fourier transform, i.e., 



m{<pf = (|-) / cos{^',<p,hka') C{h) dhk 



If this expression could be numerically computed, at least for a modest lattice 
size, we could then numerically determine \E'(0), the potential V(0), the 
lattice Hamiltonian 7i, and finally the lattice action with which to perform 
Monte Carlo calculations for the various correlation functions of interest. 

Admittedly, all that appears to be a rather formidable task. Therefore, 
one of the aims of this paper is to make this task appear more attractive. 

5 A Linear Superposition of Poisson 
Distributions 

Up to this point, we have placed all our bets on studying generalized Pois- 
son distributions as candidates for the ground-state distribution. We now 
broaden the base of our inquiry to include a much larger class of ground-state 
distributions based on linear superpositions of the very class of distributions 
we have featured so far. 

In one-dimensional probability theory, it is noteworthy that a linear su- 
perposition of even Gaussian distributions over their variance parameter is 
dense in all possible even probability distributions. This is likewise true for 
their characteristic functions; namely, one may study the set of even charac- 
teristic functions C(— s) = C{s) that may be obtained from the relation 



as / ranges over (signed) generalized functions. It is easy to see that all even 
characteristic functions and thereby all even probability distributions may 
be obtained this way.^ A similar question can be raised about linear sums of 
even Poisson distributions, or equivalently about their even characteristic 

^For example, let /(r) equal (5*^™^(t— 1) to generate the elements s^™e~^^, m = 0,1, 2..., 
take linear sums to construct the even order Hermite functions, and then note that sums 
of such functions are dense in all even functions, and hence in the subset of such functions 
that are characteristic functions. 
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functions. In equation form, this latter question reduces to the class of 
characteristic functions that can be reached by linear superposition (say) 
of the limited set of Poisson characteristic functions 



C{s) 




/(r, a;) exp{— t/[1 — cos(sii)] e du/u^} dr duj 



again as / ranges over (signed) generalized functions. Since we can recon- 
struct the Gaussian case from the Poisson case, it is clear that once again such 
linear combinations are dense in the set of all even characteristic functions 
and thereby all even probability distributions. 

In this section we study expressions for the lattice-space characteristic 
function associated with the (even) ground-state distribution in the form of 
a linear superposition of generalized Poisson characteristic functions given 



The basic ingredient in this superposition is the expression 

Cr{h) = exp{-/[l - cos(E'^/ife0jka*)]p^(0, a) n'j.d0jk} , 

which, as discussed in Sec. 4, is the characteristic function for a generalized 
Poisson distribution, an infinitely divisible characteristic function. Note well 
that we do not require that C(^K){h) is an infinitely divisible distribution itself; 
instead, we use a suitable set of infinitely divisible distributions as "building 
blocks" for C(K)ih). 

Since a limit of Poisson distributions exists that converge to a Gaussian 
distribution, we may allow for such a possibility in the continuum limit as 
a — > 0. In this sense, the expression given for our basic element incorporates 
both Poisson and Gaussian distributions, and therefore potentially includes a 
wide class of infinitely divisible distributions. With regard to the superposi- 
tion by an integral over r G MP, for some p, one may focus initially on positive 
measures K{T)dT (thus including 5-functions) such that jK{T)dT — 1. How- 
ever, it is also possible to extend such examples to suitable signed measures 



by 





drK{r) exp{-/[l 
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K{T)dT such that j\K{T)\dT < oo, jK{T)dT = 1, and > 0; we con- 

tent ourselves with a dense set of results and do not pursue examples where 
K{t) is a generalized function. Our purpose in this exercise is to increase, 
broadly, we hope, the scope of characteristic functions and thereby of ground 
state distributions that may be considered. We make no claim that our lim- 
ited family of sTipcrpositions is exhaustive of all such distributions (unlike 
the one-dimcnsional case). 

By assumption, the continuum limit of our lattice-space formulation C(^K){h) 
is intended to retain the basic structure of the separate Poisson characteristic 
functions that make up the linear combination. This implies that in the con- 
tinuum limit the resultant expression for E{K){h) is generally non-Gaussian. 
It is our proposal to look for the ground-state distribution of various scalar 
models within the family of characteristic functions offered above, and in par- 
ticular for those models that arc pscudofrcc models (thus having a limited 
number of parameters). As emphasized previously, such pseudofree models 
may serve as starting points for a perturbation analysis for more interesting 
interacting models. 

On one hand, it may even be possible to find the ground-state distribution 
for some interacting models within the family obtained by linear superpo- 
sition. On the other hand, although the class of ground-state distributions 
\['(0)^ that may be described by the linear superposition of infinitely divisible 
distributions may be relatively large, the reader may well ask by what reason 
could we expect to find an acceptable description of (say) a 99^, n > 5, model 
within that particular class of characteristic functions. 

To address this question consider the following one-dimcnsional quantum 
mechanical problem: A system with classical Hamiltonian H{p, q) is quan- 
tized by adopting the quantum Hamiltonian 

n = H{P, Q) + hY{P, Q) , 

where Y is usually present to account for possible factor ordering ambiguity. 
However, Y could in fact be a quite general operator (so long as Ti is self 
adjoint) because in the classical limit in which the original classical 
Hamiltonian i?(p, q) is recovered. This kind of ambiguity is always present 
in quantum theory, and the traditional way to deal with it is to "appeal 
to experiment". (Remark: For example, such additional terms are surely 
present in any analysis of nonrenormalizable models based on the countert- 
erms suggested by perturbation theory.) 
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It is possible that approaching a wide class of scalar models through the 
family of characteristic functions represented by C(^K){h) may correspond to 
a realization of certain models with the presence of selected counterterms 
(the analogue of Y{P, Q) above). Since a principal goal of the present study 
is to achieve nontriviality, any satisfactory set of counterterms will suffice 
initially; limitations that may arise from an appeal to experiment can be 
considered later. 

In the following section we will discuss correlation functions for the case of 
a ground-state distribution which is a generalized Poisson distribution. That 
same discussion can easily be extended to also include those distributions 
that are given by the linear superposition of Poisson distributions. 

6 Correlation Functions 

Even without fully defining our choice for a full-time, lattice-space action 
function including the required auxiliary potential, we can nevertheless draw 
some important general conclusions. In particular, let us show that the full- 
time correlation functions can be controlled by their sharp-time behavior 
along with a suitable choice of test sequences. 

6.1 Correlation function bounds for general 
distributions 

Let the notation 

denote the full-time summation over all lattice fields where {uk} denotes a 
suitable test sequence. We also separate out the temporal part of this sum 
in the manner 

Observe that the notation 0^/ (with the prime) implies a summation over 
only the spacial lattice points for a fixed (and implicit) value of the temporal 
lattice value, ko- 
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Let the notation ((■)) denote full-time averages with respect to the field 
distribution determined by the lattice action, and then let us consider full- 
time correlation functions such as 

• • • 0w(29) ) = S (1) (2) (2q) O^'^ 0u'(2) ' " " 0„'(2q) ) , ? > 1 , 

'^0 '""O '■■■'"'0 

where the expectation on the right-hand side is over products of fixed-time 
summed fields, (pu', for possibly different times, which are then summed over 
their separate times. All odd correlation functions are assumed to vanish, and 
furthermore, ( 1 ) = 1 in this normalized spacetime lattice field distribution. 
It is also clear that 

I 0u(2) • • • ) I < S (1) (2) (2q) O^'^ | 0u'(2) " " " 0u'(2g) ) | . 

'*'0 ' ' ■■■ ' 

At this point we turn our attention toward the spatial sums alone. 
We appeal to straightforward inequalities of the general form 

{ABf < {A''){B^) . 

In particular, it follows that 

(0«'(l)0„'(2)0„'(3)0„/(4))^ < (0^,(1) 0^,(2)) (0^,(3) 0^/(4)) , 

and, in turn, that 

By a similar argument, it follows that 

|(0„,a)0„,(2)---0„,(2.))|<n5l,[(05,,)]V2^ 

which has bounded any particular mixture of spatial correlation functions at 
possibly different times, by a suitable product of higher-power expectations 
each of which involves field values ranging over a spatial level, all at a single 
fixed lattice time. By time translation invariance of the various single time 
correlation functions we can assert that 



27 



which is defined at time , is actually independent of the time and, there- 
fore, the result could be calculated at any fixed time. In particular, we can 
express such correlation functions as 

Adding this relation to those established above completes a bound on the 
multi-time correlation function in terms of moments of the ground-state dis- 
tribution. 

6.2 Correlation function bounds for Poisson 
distributions 

Observe that the bounds we have discussed up to this point apply to arbi- 
trary ground-state distributions Let us now specialize to generalized 
Poisson distributions with which we shall be able to make additional and 
more specific remarks. 

For Poisson distributions, we note that the moments of the ground-state 
distribution for the correct (but unknown) ground-state may be directly re- 
lated to the truncated moments for the ground-state distribution, which in 
turn are directly determined by the "moments" of the weight function a). 
Let us recall the notation introduced earlier, namely 

which are not moments but noments (since J p{(f), a) H^ii^fc = oo). These no- 
ment expressions coincide exactly with the truncated moments of the ground 
state distribution. Recall that ordinary and truncated moments are related 
by the generating function relation 

(e""^"') = exp[(e°'^"' - 1)^] , 

which holds for all a, where T denotes truncated. Therefore, in our case. 

For example, we list a few relations that follow from this connection: 
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(0t')+3(0^o^ 

(0^0 + 15(0^0(0^) + 15(0^0', 
(0^') + 28(0^,)(0!')+35(0^,r 



+ 21O(0^,)'(0t') + lO5(0^,)\ 



etc. 



The conclusion of this exercise is that all of the multi-time correlation 
functions are bounded by terms that are composed from the noments deter- 
mined from the weight function p{(f),a). If we can choose p{(f),a) so that it 
has a suitable continuum limit, and in such a way that the noments remain 
finite in that limit, then we will have obtained a bound on every full-time 
correlation function in the continuum limit. Convergence (of subsequences, 
if necessary) of these correlation functions follows. Moreover, the ground- 
state distribution stays firmly within the family of Poisson distributions and, 
by working with the various noments (0^?), we have avoided the fate of a 
Gaussian, i.e., trivial, limiting distribution. 

7 From Noments to Genuine Averages 

7.1 A two-point normalization 

The discussion regarding correlation functions has shown that many func- 
tions of interest can be evaluated exactly, or are suitably bounded, by terms 
involving various noments such as 



where q > 1 and Uk denotes a test sequence here used in space alone. Ex- 
pansion of the power 2q leads to expressions of the sort 



It is these noments that we wish to discuss in this section with the possible 
aim of approximately evaluating them using Monte Carlo techniques. 
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In order to invoke Monte Carlo calculational methods, however, it is clear 
that we need to introduce a normalized probability measure, and for this 
purpose it will suffice to fix the factor R{a), at least implicitly. We are free 
to fix this normalization based on the freedom involved in the usual field- 
strength renormalization. And we do so by selecting a carefully chosen and 
suitably weighted two-point function and then declare that the noment of 
this expression is unity. 

To this end, let us choose Yk^i - where, in this case, k and I involve only 
the spatial components - as the elements of a symmetric matrix Y = {1^,/} 
which is positive definite. Furthermore the elements Y^,/ are chosen so that 
their "high frequency" components are suitably reduced. We will offer two 
examples of what we deem to be suitable matrices, and from those it should 
become clear just what we have in mind regarding the "high frequency" 
components. The first example we choose for ^ is suitable for a lattice 
formulation of the problem, such as may arise in Monte Carlo calculations, 
where the lattice size L and the lattice spacing a may vary, but remain 
finite, and no attempt at approaching the continuum limit is contemplated; 
this example will be used in Sec. 7.1. The second example we choose for Ykj 
is one designed to be use for an approach to the continuum limit; the second 
example will be used in Sec 7.2 in which the continuum limit is discussed. 

Observe that the unit matrix 6k,i is positive definite, but it does not 
discriminate between "low" and "high" frequencies. To describe the unit 
matrix in other terms, we first assume that L, the number of lattice points 
along each axis, is odd, or in other words, L — 2P -\- 1, where P is a positive 
integer. We choose two integers r and s, where —P < r < P, and similarly 
for s, and observe that 

^ g27ri(r-s)m/L = ^ _ 
^ m=-P 

Now, as a preliminary to making our choice for Y^^i, we dampen the weight 
of the high frequencies in the previous sum by defining 

Z^^=— ^ ^-\m\T ^2-Ki{r-s}m/L ^ 
^ m=-P 

where T > is a damping parameter to be chosen. Clearly, the resultant 
function is still positive definite as desired. Indeed, the sum can be evaluated 
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in closed form as 

{1 - - 2e-(^+i)^ cos[27r A (P + 1) / L] + 2e-(^+2)^ cos[27r AP/L]} 
L{1 - 2e-^ cos[27rA/L] + e'^^} ' 

where A = r — s. Finally, to define our first cfioice for Y^^;, we let 

s 

i=i 

where, as noted above, k and / involve only the spatial components. If, on 
the other hand, L is even, i.e., L = 2P, for integral P > 0, then we choose 

Z^^ = — ^ g-|m|Tg27ri(r-s)m/L ^ 
^ m=-P+l 

which has the closed form 

7 _M .-PT, ,^A1r 2[l-e-^cos(27rA/L)] 

^,,,-[1 ^ U[l-2e-^cos(27rA/L) + e-2r] J' 

The quantity Y^,/ is constructed from Z^^s in the same way as for the odd L 
case. 

Although not needed immediately, we also choose our second form for Yk^i 
suitable for discussing the continuum limit. For the following, let z G M, and 
let the ortho normal set of Hermite functions be called hn{z), < n < oo, 
with the usual property that 

oo 

J2hn{z")hn{z')=6{z"-Z'). 
n=0 

We dampen the high frequencies once again and define 



Z{z",z') = Y,^-''^hn{z")K{z'), 



n=0 

where again T > is a parameter to be chosen. This sum can be evaluated 
in closed form as 

Z{z", z') = -L= exp[-(l/2)(/'2 + z'^) coth(r) + z"^'csch(r)] . 
V 27rT 
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Clearly, the kernel Z(^z" ^ z') is positive definite. Finally, we define the function 
Y{x, y) - where both x = {xi, X2, ■ ■ ■ , Xs} and y = {yi, 1/2, ... , ys} are in 
- by the expression 

s 

Yix,y) = YlZ{xj,yj) . 

i=i 

Note that since each Hermite function qualifies as a test function, we may 
refer to Z{z",z') and Y{x,y) as positive definite "test kernels". The given 
choice Y{x, y) is suitable for the continuum limit for the whole space M*, but 
if one chooses to use it for a large lattice, before taking the continuum limit, 
it would be acceptable to use 

Yk^i = Y{ka, la) , 

namely, the values that Y{x,y) assumes on the lattice points themselves. 

Having chosen a suitable, positive definite matrix as discussed above, 
we are in a position to establish our normalization criterion. In particular, 
for a given, fixed, choice of the matrix Y, we (arbitrarily) declare that the 
noment 

This normalization condition fixes the constant R{a). In effect, we are saying 
that a certain two-point function with a suitably chosen positive definite 
weighting has a noment of unity. We can put this normalization to good use 
for us as follows. 

For all > 1, we rewrite the noment of general interest as 



J n',[E;j,_i<A2]i/2 '-V 
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This last equation establishes the noment of interest as a genuine average 
involving a normalized probability measure 



da{(j)) 



[E'(f)kYk,i(l)ia^']e 



X 



n' 



3fe • 



It is the final expression for (0;^ 4>i2 ■ ■ ■ 4'hq ) "^^ would like to evaluate approx- 
imately by means of conventional Monte Carlo methods. 

We have arrived at the expression of noments as genuine averages on 
the basis of a convenient choice of normalization given above. Note that if 
instead we had chosen to normalize matters so that 



Ria) 



rather than 6 = 1 as above, then the final result would have been 

{<t>h<t>h ■ ■ ■ <t>h,) [ 't\ ^^('^) ' 

where the normalized probability measure c?(t(0) is unchanged. Hereafter, 
we assume that 6=1. 

It should be noted that there is a special advantage derived from the fact 
that the noments can be expressed as suitable averages over a particular prob- 
ability distribution. This advantage follows from the fact that inequalities 
such as 



/ 



{<t>h <t>h) 



2s 



da{(p) > 



{<t>h <t>h) 



T,'4>kYn,i<t>i(^''' 



da{(p) 



hold. 



7.2 Continuum limit 

Finally, we briefly take up the question regarding the continuum limit. There 
are several important and distinct aspects of the continuum limit, namely. 
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not only a but N' ^ oo, and in fact, it is also necessary that the spacial 
lattice volume iV'a* — > oo. Let us first discuss the situation for large N'. 

When dealing with a large number of integration variables, i.e., N' ^ 1, 
it is important to note that there are important - even profound - differences 
between Gaussian-like integrals and Poisson-like integrals. In presenting the 
following discussion wc closely follow a section of [3] dealing with an idealized 
set of examples. For all p > 1, consider the two sets of integrals 

as representative of the Gaussian-like expressions, and 

Ip{2p) = J {K4re-^'''>'^HK4]-'''^'KdA , 

as representative of the Poisson-like expressions. These integrals arc "cari- 
catures" of the ones we study, but the results are nevertheless informative. 

To help in our study, let us introduce hyper- spherical coordinates [3] (c.f., 
also [9]) defined by 

(j)k = Krjk, < K < oo , -1 < r^fc < 1 , 

Here k denotes an overall radius variable while {rjk} denotes an iV'- dimensional 
direction field. In terms of these variables, the Gaussian-like integrals take 
on the form 

Ig{2p) = 2 J K'^e-^^\'''-' dKS{l - J:Wk) KdVk ■ 

Evaluation of this expression for large A^' is dominated by the factor k^'~^, 
and a steepest descent method can be used to evaluate the integral over k 
to a suitable accuracy. To leading order, it follows that the stationary point 
is given by «; = [N' /2AYI'^. As a consequence, for each value of A, the 
integrand is supported on a disjoint set of k as A^' — oo. This well-known 
fact leads to divergences in perturbation calculations. For example, let us 
calculate 
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for a different value of A by the perturbation series, 

7^(2) = 7g(2) - AA 7g(4) + ^AAf 7^(6) - • • • , 

where AA = A* - A. Since 7g(2p)/7g(2) oc N'^p-^\ this series exhibits 
divergences as A?"' — > oo. 

Let us now consider the Poisson-hke integrals expressed in the same 
hyper-spherical coordinates. It follows that 

7p(2p) = 2 J K^Pe-^^' [«:2]-^72/,iv'-i ^^§{1 - U'drjk ■ 

Here we see no large K-power in the integrand, and it follows, e.g., that 

7p(4) _ jK^e-^'^' dK _ 1 
7p(2) ~ /Ace-^«' dK ~ A ' 

which is finite and independent of N' . Hence a perturbation calculation of 
the change of A to A* for Ip{2), for example, involves no divergences as 
A^' ^ cxd! Further examples of the difference of Gaussian- and Poisson-like 
integrals is offered in [3], but the main point has already been made showing 
the profound difference between these two types of integral for large N'. 

How do these examples impact on our main discussion? Consider the 
expression for a simple noment given by 

Expressed in hyper-spherical coordinates, this integral becomes 

A: L I k ^ ll \ 

and we see an analogous absence of the parameter N' regarding the integral 
over K. This favorable situation has arisen from the presence of a denominator 
factor which, although different in detail, has an overall power of k identical 
to that of the elementary examples for Ip{2). Ideally, one might like to 
choose a denominator factor that is as local as possible on the lattice, and 
that would mean choosing a set of {Pk} parameters of the form (3k = P 5fc,o- 
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This choice has the desired locahty property, but it does not lead to an 
integrable denominator, as noted previously and which is further discussed 
in the Appendix. The choice of Pk — Jk - that is, our featured fcP'"^ set of 
values - is a suitable compromise that has locality in the continuum limit but 
retains enough lattice spread to always ensure an integrable denominator. In 
calculations of our principal interest, and with our featured choice of the 
denominator factor, we expect to benefit from the same lack of divergences 
due to the missing large TV' dependence of the hyper-spherical radius k as 
was exhibited in the elementary examples. 

Next we consider the problems associated with the limit a ^ along 
with N' ^ oo, ideally, so that N'a'"^ oo, but at the very least, so that 
N'a'^ — > V, where < V < oo, and preferably V is "large" in some sense. 
To appreciate the kind of problems faced it is useful to examine first the 
usual situation for a "free" theory and its continuum limit. 

Consider the lattice formulation for the ground-state distribution of the 
free (= Gaussian) theory given by 

where M denotes a normahzation so that Cf{0) — 1. It is clear in this case 
that the two-point function is given by 

i.e., essentially the r, s matrix element of B~^, the inverse matrix to B. A 
suitable continuum limit of this expression starts by considering the smeared 
expression 

for sequences {^r}, for example, such that Tl^u^a^ — 1. In this case, the 
continuum limit should assume the form 

{{ju{x)(t){x) dfxf) = \ju{x)u{y)F{x,y) d^d\i , 

for some generalized function F{x,y) and for all functions u{x) G C^, for 
example, for which Ju{x)'^ d^x = 1. The key concept is the limiting behavior 
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of the lattice two-point function to a suitable generalized function in the 
continuum limit, a well-known property for a free theory. Observe that the 
properties of F{x, y) are directly determined by the properties of the matrix 
B~^, the inverse of the matrix B. 

Now let us consider a similar limit for the ground-state distribution de- 
termined by a generalized Poisson distribution. In particular, we first focus 
on the two-point function 

where, as discussed above, normalization is set by 

As in the Gaussian example, we now consider the expression 

and, for a satisfactory continuum limit, we ask that as a ^ and N'a^ oo, 
or at least N'a^ V 1, we find that 

{{Ju{x)(l){x) d^xY) = Ju{x)u{y)G{x,y) d^xd^y , 

for a suitable generalized function G{x,y), say, for the same set of smooth 
functions u considered above. In particular, according to our normalization 
procedure for determining R{a) , we have required that the continuum choice 
of the positive definite test kernel, Y{x,y), such as the example introduced 
in Sec. 7.1, satisfy 

jY{x,y)Gix,y)d'xd'y^l. 

Note that, quite unlike the Gaussian case, the matrix A is only indirectly 
involved in the determination of the generalized function G{x, y). Indeed, to 
have a suitable form for G{x, y) , there is no guarantee that the matrix A has 
any clear continuum limit by itself, i.e., there is no requirement that 

T.'^^i(t)kAk,i(t)ia^' j(t){x)A{x,y)(t){y) d'xd'y 
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for a suitable generalized function A{x, y) in the continuum limit. All that 
is required is that G{x, y) be suitable - as well as suitable distributional 
behavior for all the higher-order correlation functions as well! 

It would be helpful to build up some computational experience for correla- 
tion functions for the generalized Poisson distributions in order to eventually 
help select matrices such as A by connecting them to the generalized (corre- 
lation) functions (such as G) to which they give rise. Perhaps Monte Carlo 
calculations could make a contribution to that effort. 

In summary, in our all-too-brief discussion of the continuum limit, we have 
observed: (i) the atypical behavior of the limit in which N' ^ oo, namely, 
in the case of generalized Poisson distributions, how, in hyper-spherical co- 
ordinates, the appearance of A^' as a large power of the overall field radius is 
absent, thereby rendering perturbation series termwise finite in contrast to a 
comparable calculation for Gaussian-like integrals; and (ii) the different role 
of internal parameters (e.g., matrices A and B) in the ground-state distri- 
bution for Poisson- and Gaussian-like distributions, and how they effect the 
form taken by the continuum limit of correlation functions. 

8 Conclusions 

In discussing lattice quantum field theory models, it is conventional - and, of 
course, very natural - to start with a choice for the lattice action. A lattice 
action for a model such as Lpf^, n > 5, however, necessitates that some difficult 
choices must be made for the appropriate counterterms. Appropriate in this 
sense means counterterms that have the virtue that the resultant model and 
its continuum limit will faithfully describe the original model chosen, as well 
as contain the original motivating classical model in the limit that ^ — >■ 0. 
This is a lot to ask, and it is not surprising that the auxiliary potential (i.e., 
counterterm) generally chosen for this purpose does not perform as desired. 

In this paper we have taken a different approach. First of all, recogniz- 
ing that counterterms suggested by regularized, renormalized perturbation 
theory are inappropriate for nonrenormalizable models (see Sec. 2). we have 
felt compelled to direct our initial focus away from the lattice action and 
place it instead on the sharp-time, ground-state distribution function ^'(0)^. 
To avoid triviality, this distribution must be non-Gaussian, and to avoid the 
basin of attraction that leads to a Gaussian distribution (and its associated 
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triviality), we have chosen to emphasize the possible relevance of the only 
other class of suitable infinitely divisible distributions, namely, the gener- 
alized Poisson distributions. Such distributions avoid manifest triviality by 
being non-Gaussian, and they are robust to survival - as are the Gaussian 
distributions - under general limiting conditions. 

We found that ground-state distributions in general, and Poisson dis- 
tributions in particular, were faithful replacements for the lattice action 
which implicitly - even though somewhat indirectly - was determined by 
the ground-state distribution. Moreover, and this is the important point, 
the class of lattice actions associated with generalized Poisson distributions 
are automatically assured to lead to manifest nontriviality by ensuring a 
non- Gaussian continuum limit. 

Although we could be fairly specific about which Poisson distribution 
we favored, it was nevertheless extremely difficulty to see just which lattice 
action was implicitly connected with that distribution. Thus we felt necessary 
to abandon any focus on individual models, and rather focus on classes of 
models hoping that such a wider net might capture one or another model of 
possible interest. 

In the putative continuum analysis of Sec. 3.4, it was emphasized that 

suitable matrix elements of the canonical momenta, as well as those of the 
Hamiltonian, were uniquely specified by our choice of the characteristic func- 
tion associated with the ground-state distribution. It is interesting to further 
note that if it were possible to actually determine the associated canonical 
coherent states defined with the ground state itself serving as the fiducial 
vector, then additional matrix elements of interest could be computed. In 
particular, incorporating and extending the notation for \h) of Sec. 3.4, we 
let 

\h,g) = e'iJ[K^)4'{^)-9{x)*i^)]d''^ |o) 
for which it follows [10] that 

{h, g\4>(x)\h, g) = g{x) , (/i, g\T^{x)\h, g) = h{x) , 
and more significantly (assuming a model under consideration) that 
hm (/i, g\ H \K g) = /{ \h{xf + \[Vg{x)f + Im-" g{xf + \g{xY}d^x , 
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implying the close connection of the diagonal coherent state matrix elements 
of the quantum Hamiltonian with the classical Hamiltonian itself [10]. All in 
all, these conditions place strong limits on the ground-state wave function, 
even if they are highly implicit. 
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Appendix: Fundamental Restrictions on the 
Choice of the Weight Function a) 

In this Appendix we address the question of acceptable values for 7 and 
acceptable choices for the set {(3k} that are both part of the canonical form 
for a) we have adopted in Sec. 4. Therefore, we again consider the 
integral given by 



which should diverge when p = and converge whenever p > 1. We as- 
sume that f/(0, a) controls large {(pk} behavior for all p, and so the question 
of divergence or convergence is related to the small {(f)k} behavior where 
U{((),a) ~ 0, and where U{4>,a) is not a contributing factor. Indeed, for the 
general question of divergence or convergence the specific form of U (0, a) is 
not of particular relevance and it may effectively be replaced by 




U{^,a) = E',^la\ 



leading to the expression of interest 
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It is convenient at this point to again introduce hyper-spherical coordi- 
nates for which 0^ = Krjk, n > 0, —1 < rjk < I, where SJ,?]^ = 1, and 
^'k^k — i^"^- III terms of these coordinates, the integral of interest reads 



Ip^2R{a) J 



For p = we attribute divergence of this integral to a divergence at k = 0. 
Consequently, we need 2^N' > N', i.e., 7 > 1/2. For p > 1, convergence at 
K = requires that 2-fN' < N' + 2, i.e., 7 < 1/2 1/N'. Thus we are led to 
the bounds 

1 1 1 

- < 7 < - H . 

2 - ' 2 N' 

Since the limit iV' ^ 00 is eventually to be taken, we already satisfy this 
situation by adopting 7 = 1/2, as noted earlier. However, we are not quite 
finished because convergence for p > 1 also requires convergence of the inte- 
grals over the {rjk}, and as we now shall see, this convergence requires some 
restriction on the set of coefficients {Pk}- 

The first remark regarding such integrals is to observe that with 7 = 1/2 
it is not possible that f3k oc 6ko, for in that case integrals over each rjk, save 
when k = t, lead to a divergence. Consequently, Pk must be nonvanishing 
for additional sites. To study how many additional sites are sufficient, it is 
convenient to express our basic integral in yet another form, namely 



= R{a) a^^'/2^-iv'/2 J J 0^2P[n;,^(-V2)] e-s;(i+s;/3.-<A,)<A^^ U',dXkKdcl>k , 

where the integral for each variable runs from to cxo. Exchanging the 
order of integration leads to 



Ip = K 
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where 



K = 2-P{2p-l)\\R{a) a 



—sp 



Observe that when the basic integral Ip is expressed in the {Xk} variables, 
divergence or convergence concerns the behavior of the integrand for large 



Initially, let us again show, when p = 0, that I — oo - this time in terms 
of the {Xk} integration variables. This result follows from the lower bound 
given by 



which holds because f3k < 1- If we introduce coordinates of the form = ^77^, 
C >. 0, —1 < rjk < 1, and Sj,7y^ — 1, then it follows that 



The first integral diverges, while the remaining factors are positive. Thus we 

have again shown that /q = 00, as expected. 

For p > 1, wc require that Ip < 00. First wc observe by construction that 
Iq < Ip whenever q > p. Therefore, to show that /p < cxd for p > 1, it suffices 
to show that /i < 00. Thus we focus attention on 



If this integral were to diverge it must diverge for large A. We can reach 
large A in many different directions. In particular, we can imagine that P, 
1 < P < N', oi the A variables are all approaching infinity in a certain 
direction. As an example, consider P — 3 and Ai = .5^, A2 — .2(, X3 — .1^, 
as C — > 00, while all other A variables have finite values. The "direction" 
of approach to infinity in this case refers to the relative size of the growing 
A terms, i.e., the direction as determined by rjf = .5, r/l = .2, and r/l = .1. 
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Staying with this example for a moment, the three large A values will arise 
in several terms in the denominator because of the possible many-fingered 
nature of the /3 terms. We divide the set of distinct space-like k values into 
the set P which contains one or more of the P variables approaching infinity 
and the set Q containing the remaining N' — P terms for which the variables 
A do not approach infinity. In dividing the appropriate sets below, we have 
added an asterisk to the set of A^ values that are approaching infinity to 
serve as a reminder of that fact. Thus we are led to the expression 



In the numerator of this expression there are P variables approaching 
infinity at the same time, where 1 < P < N', and we need to consider all 
possibilities. For the denominator, there are several cases to consider. We 
assume that the set {Pk} has 1 -|- 5*0 nonvanishing terms and that they are 
connected to the "home coordinate" , k — 0. In particular, let us focus on the 
choice fcP'"^ (consisting of the point k and its 2s nearest neighbors in spatial 
directions) for which 5*0 = 2s. As a consequence, there are large A variables 
in P -|- (+1) factors in the denominator; here the term S, which may vary 
from case to case, arises from the many-fingered nature of {Pk}, and the 
term (-1-1) adds 1 only if a large A factor appears in the denominator factor 
associated with p = 1. Convergence of Ji is ensured if for all P, 1 < P < iV', 
the choice of {f3k} ensures that S'(+l) > 1 in all possible combinations of 
the large A values. A few specific examples may help clarify how the choice 
of /cP'"^ fulfills all the necessary requirements. 

For example, if P = 1 then S (+1) = Sq (+1) > 1. If P = A^' - 1 includes 
all the lattice points except one, then 5'(+l) = 1; if P = A^', and thus 
includes all lattice points, then 5'(+l) = 1. If P = L, all of which lie in a 
single coordinate direction (with period L), then S (+1) = 2(s— 1)L (-1-1) > 1. 
Incidentally, this example shows that fingering in just one spatial direction 
is not sufficient - and, by analogy, fingering in (s — 1) spatial directions is 
also insufficient. It is necessary to finger in all s different spatial directions. 

A picture may also help; see Fig. 1. Although not directly relevant to 
nonrenormalizable fields, we give a two-dimensional example (s = 2) for 




X 



[n',(i + E;,p/3fc_Ar + s' A_,A0V2] ■ 
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which the points with large A values arc marked with solid circles •, the 
additional points in which those large values appear in denominator factors 
are marked with an open circle O, and the point k — t is marked with an x. 
In the example of this picture, P = 2, Sq = 4, S = 6, and since the p — 1 
term for this example does not contribute, the (+1) term is absent. 
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Fig. 1: A two-dimensional example with two large A values (dark circles) 
which also have a presence in six additional locations (open circles) . The 
place marked by x is where the additional moment is located. 
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